Approach to Universality in Axisymmetric Bubble Pinch-Off 
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The pinch-off of an axisymmetric air bubble surrounded by an inviscid fluid is compared in four 
physical realizations: (i) cavity collapse in the wake of an impacting disc, (ii) gas bubbles injected 
through a small orifice, (iii) bubble rupture in a straining flow, and (iv) a bubble with an initially 
necked shape. Our boundary-integral simulations suggest that all systems eventually follow the 
universal behavior characterized by slowly varying exponents predicted in [Eggers et al, PRL 98, 
094502 (2007)]. However, the time scale for the onset of this final regime is found to vary by orders 
of magnitude depending on the system in question: while for the impacting disc it is well in the 
millisecond range, for the gas injection needle universal behavior sets in only a few microseconds 
before pinch-off. These findings reconcile the different views expressed in recent literature about the 
universal nature of bubble pinch-off. 

PACS numbers: 47.55.df, 47.15. km, 47.11. Hj 
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The precise nature of axisymmetric bubble collapse in 
a low- viscosity fluid has been a subject of controversy 
over the last years. Such a collapse may be initiated 
by a variety of different forces (e.g. surface tension, hy- 
drostatic pressure, and external flows). In a later stage, 
however, it is only the requirement of mass conservation 
that forces the liquid to accelerate more and more as 
the shrinking bubble neck closes in on the axis of sym- 
metry. This purely inertial nature of the final collapse 
motivated the first hypotheses about universality of the 
final collapse regime [l], Q . A simple power-law was pre- 
dicted with the neck radius scaling as the square root of 
the time remaining until the pinch-off singularity. Nei- 
ther numerically nor experimentally could this behavior 
be confirmed. Instead, for different systems and initial 
conditions a variety of scaling exponents all sli ght ly above 
1/2 have been obtained [SEEESEEESEll] leading 
to doubts about the universal nature of bubble collapse. 

Recently, the idea of universality has been revived by 
[l2l . EH who suggested an intricate coupling between the 
radial and axial length scale. The authors of [TH ex- 
plicitly predict the existence of a final universal regime 
which however is no longer a simple power-law, but char- 
acterized by a local exponent that slowly varies in time. 
The value of 1/2 is recovered in the asymptotic limit in- 
finitesimally close to pinch-off. According to this theory 
the variety of observed exponents corresponds to different 
time averages of this local exponent. 

In the present work we aim to reconcile the different 
views about universality in axisymmetric bubble pinch- 
off expressed over the last years. The key aspect is that 
we examine in detail how and when different physical re- 
alizations of bubble pinch-off reach the universal regime. 
We present detailed numerical simulations which are able 
to follow the neck evolution over more than 12 decades 
in time even for complex realistic systems. We first show 
that all systems that have recently been studied in the 
context of bubble pinch-off eventually follow the same 
universal behavior predicted by jl3| . The time scale on 



which universality is reached, however, varies enormously 
from one setup to another. For an impacting disc [S E3| 
universality can be observed during several milliseconds 
prior to pinch-off and thus on a time scale which is exper- 
imentally accessible. However, for gas bubbles injected 
through a small needle 0, E i, E 0, E E ES EH, EEL El 
universality sets in only a few microseconds (or even less, 
depending on the precise initial conditions) before pinch- 
off. This may well be the reason why universality has 
thus far never been observed even in very precise gas in- 
jection experiments and why non-inertial effects like sur- 
face tension have been claimed to play a dominant role 
in this geometry [IS EI| • 

Four different physical systems have been reported in 
the literature on bubble pinch-off, numerically and ex- 
perimentally, and will be compared in this Letter: 

(i) Impacting disc. The bubble is created by the impact 
of a circular disc on a liquid surface E EH as shown in 
Fig. []Ji). Upon impact an axisymmetric air cavity forms 
and eventually pinches off halfway down the cavity under 
the influence of hydrostatic pressure. Immediately after 
pinch-off the ejection of a violent jet can be observed 
whose formation however is not caused by the singularity 
alone [l7j as one might naively expect. Since surface 
tension is negligible E E3, E3 the only relevant control 
parameter is the Froude number Fr = Vq /gRo with the 
impact velocity Vo , gravity g and the disc radius Rq . In 
the data reported here the disc radius varies between 1 
and 3cm and the impact velocity ranges from 1 to 20m/s. 

(ii) Gas injection through an orifice. A small nee- 
dle sticks through the bottom of a quiescent liquid pool 
[E E 1 E I, E ES, EH, El, UM as illustrated in 
Fig. [TJii) . A pressure reservoir connected to the needle 
slowly pushes a gas bubble out of the needle's orifice. The 
bubble then rises under the influence of buoyancy. When 
the air thread between the orifice and the main bubble 
becomes long enough, surface tension causes the thin- 
ning of the neck which eventually leads to the pinch-off 
of the bubble. We present data for three sample config- 
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FIG. 1: Illustration of the bubble collapse in the four dif- 
ferent systems: impacting disc (i), gas injection through a 
needle orifice (ii), bubble in a straining flow (iii), and initially 
necked bubble (iv). Blue lines correspond to the free surface 
at pinch-off, while dashed and dotted black lines represent 
earlier bubble shapes. The disc and the needle are depicted 
in red. 



urations A-C corresponding to Figs. 4, 10, and 6 of 
and characterized by a Weber number WeA,B,c =0.007, 
36 and 173, respectively. (The Weber number is defined 
as We = pQ 2 / (ir 2 a 3 a) with water density p, gas flow rate 
Q, needle radius a and surface tension a). 

(iii) Bubble in a straining flow. The initially spherical 
bubble collapses due to a surrounding hyperbolic strain- 
ing flow @, G1 [HIS, see Fig-ffiiii)- 

(iv) Initially necked bubble. Surface tension causes the 
pinch-off of a bubble starting off with an initially already 
pronounced neck [l3[ as illustrated in Fig. QJiv). 

In all four systems we consider the idealized invisicid, 
axisymmetric bubble pinch-off ne glecting the influence of 
the inner gas dynamics @, ©, QIH El M, HI, 



0, fiol . EH , ano - non-axisymmetric perturbations 
16J. For our numerical investigations we employ an 
axisymmetric boundary-integral code similar to the one 
described in [l4| which has shown very good agreement 
with experiments of system (i) for various impact geome- 
tries 

HQ HI- The validit y of 

our implementation for 
the other systems is verified by comparison with the bub- 
ble shapes obtained in various earlier works @, EH, EH • 
Some details about the simulation parameters are given 
in the EPAPS file [H]. 

In a first approach to an analytical description of bub- 
ble collapse, the bubble shape can be approximated as 
an infinitely long cylinder (neglecting axial velocities) 
which yields a two-dimensional version of the well-known 
Raylcigh equation 0, [H, H, [H, E3 for the neck radius tq 



Here, F represents the force driving the collapse, and 
dots denote the derivative with respect to time t. R^ is 
a cut-off radius required to saturate the pressure at large 
distances. Assuming R^ to be a constant leads to the 
neck radius r$ shrinking as a power-law (possibly with 
logarithmic corrections [3. lol Il2| ) with exponent 1/2 as 
a function of the time to pinch-off r = t c — t. where t c is 
the closure time. At first sight, this expectation seems to 
be very well confirmed for all four systems by the lines in 
Fig[21 (b) which to the naked eye appear perfectly straight 
over more than 12 decades. The slope which corresponds 
to the scaling exponent is slightly larger than 1/2, in 
agreement with previous experiments and simulations [H, 

5HEi,0,BQQ,El. 

A more detailed look at the local exponent, defined as 
the slope in Fig. 12(b), a(r) = 91nr /51nr (Fig. [2(c) - 
(e)) however, reveals that the behavior of the neck radius 
cannot be described by a simple power-law. The local ex- 
ponent a varies during the approach to pinch-off [l3| . In 
fact, the relevant equation for the time-dependence of a 
in [13( can be derived directly from Eq. by letting 
Roo = 2y / 7o7v Here r c is the local axial radius of cur- 
vature (see Fig. 2(a)). The combination ^fr^F c is the 
scale by which the axial coordinate has to be rescalcd 
in order to collapse neck profiles at different times when 
rescaling radial coordinates by r§ 0, 0]. This leads to 
the aspect ratio of the cavity naturally being defined as 
7 = To/ y/r r c . With the above substitutions and working 
in the limit of vanishing F we obtain from Eq. (fTj): 



da 
dlnr 



2a' 



(2) 



d (rprp) , r 
dt Roo 



F. 



(1) 



which is exactly identical to Eq. (4) in [13[ (being Fi = 8 
[23l ] and <2q = 2j 2 in the original notation). Equation 
© with the da Id In r term neglected due to the slow 
variation of a [26j is shown as the dashed red line in 
Figs. [2] (c)-(e). It represents the universal regime where 
the only driving is provided by inertia and all external 
forces have become negligible. We will now proceed to 
compare the approach of the different systems (i)-(iv) to 
this universal curve. For this we use the aspect ratio 7 
as a universal "clock" replacing the time to pinch-off r 
[23l ] . Note that 7 — > as r — > meaning that the cavity 
becomes more and more slender @, E3 • 

One of the key points to address is if and how this be- 
havior can be observed experimentally. Besides the obvi- 
ous difficulty of obtaining a sufficient number of decades 
to observe the slow variation of the local exponent, the 
crucial question is: at what time (before pinch-off) does 
the system exhibit universal behavior? This is crucial be- 
cause, firstly, the duration of the universal regime needs 
to be within the time resolution of the experimental 
equipment. And secondly, the onset of universality needs 
to happen before other effects such as air flow, viscosity, 
non-axisymmetric instabilities, etc. unavoidably destroy 
the purely inertial regime. We will now provide those 
time scales for the various systems based on numerical 
boundary-integral simulations which do not have these 
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FIG. 2: (a) Illustration of the cavity surface (blue line), the minimal neck radius rg and the local radius of curvature r c . (b) 
"Classical" plot of the neck radius versus the time to pinch-off for system (i) with Fr=5.1 (blue), system (ii) in setup B (red), 
system (iii) (dark gray) and system (iv) (magenta) . The curves are offset by one decade from each other to improve readability, 
(c) The local exponent a as a function of the aspect ratio 7 for system (i) with Fr=3.4 (cyan), Fr=5.1 (blue), Fr=46 (brown), 
Fr=500 (green), and Fr=4000 (black). After an initial transient all curves follow the same universal regime. The dashed line is 
Eq. (J2J). The local maxima correspond roughly to the start of the universal regime, (d) The local exponent for system (ii) in 
the three configurations: A (red), B (gray) and C (yellow). All curves A-C lie practically on top of each other, (e) The local 
exponent for system (iii) in dark gray and system (iv) in magenta follows the same universal behavior close to pinch-off. Small 
jumps in the data are due to the crossover between different node positioning algorithms employed in the initial and the final 
stages of the simulation (see EPAPS document [25]). 



limitations. 

We start by considering the impacting disc system (i) 
in Fig. [2] (c). It is evident that the data for all values of 
the control parameter follow - after some initial transient 
- the same universal curve in excellent agreement with 
Eq. ([2]). Figure [2] (c) further gives us a good measure at 
what aspect ratio the universal regime is attained: ap- 
proximately after passing their respective local maxima 
all curves follow the same behavior. The aspect ratio of 
this maximum can then easily be related to the phys- 
ical time before pinch-off t u . We find t u m 6ms and 
t u w lms for Fr=3.4 and Fr=4000, respectively. That 
the high Froude case reaches universality later can be 
intuitively understood: at high Froude the cavity closes 
deeper and therefore the hydrostatic driving pressure is 
larger and its effects on the neck dynamics can be felt 
longer. It is remarkable nevertheless that the duration of 
the universal regime changes only by a factor of less than 
10 while the corresponding control parameter varies over 
three orders of magnitude. At the same time both values 
are easily within experimentally accessible time scales. 

We now compare this to system (ii) jtlie bubble injec- 
tion through a small needle 0, i, Hi, 0, S S M, EH, M 
d] in Fig. [U (d). The approach to the universal regime 
is much less abrupt than in system (i). Due to this more 
gradual approach, it is difficult to specify precisely the 
time when universality is reached for the gas injection 
needle. We thus choose to keep our previous definition 
of t u being the time corresponding to the local maxi- 



mum in Fig. [5] (d). This gives a good upper bound for 
the time when universality sets in. Surprisingly, we find 
even these times to be of the order of 5/zs in case A, 60ns 
in case C and as low as 10ns in case B, respectively p7j . 
Thus, the duration of the universal regime in the needle 
setup is dramatically (by at least three orders of magni- 
tude) shorter than for the impacting disc. This may well 
explain why, among other experimental limitations, the 
universal regime has never been observed experimentally 
for this widely used system. 

Figure [2] (e) confirms that also the systems (iii) and 
(iv) follow the universal regime. System (iii) does so 
even over the entire plotted range. Both are somewhat 
idealized systems and we are not aware of any experimen- 
tal investigations regarding the approach to pinch-off in 
these systems. Without a relevant length and time scale 
it is therefore impossible to specify the physical time to 
universality in these cases. 

The different behaviors of the individual systems can 
intuitively be understood as follows. System (iii) con- 
tains no external driving force other than liquid inertia 
which makes it the ideal system to compare with Eq. |2} . 
Indeed, this entirely inertial system follows the univer- 
sal regime over the widest range in aspect ratios of all 
systems studied. Similarly, the very violent disc impact 
should introduce a large amount of inertia in system (i) 
which consequently follows the universal regime also for 
a rather long time. On the other hand, the two systems 
governed by surface tension (ii) and (iv) contain little 
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FIG. 3: (a) The local Froude number for system (i) and an im- 
pact Froude number Fr=3.4 (cyan) and Fr=4000 (black) and 
the local Weber number for system (ii) with configuration A 
(red) and B (gray). The onset of the universal regime can be 
located roughly after the respective non-dimensional quanti- 
ties have become larger than order unity (horizontal dashed 
line), (b) The averaged exponent measured in [|| (black di- 
amonds) is reproduced very well by our model with the con- 
stant C — 0.7 (red line). 



inertia and thus approach the universal regime only rel- 
atively late and in a similar fashion. 

To make the above arguments more quantitative we 
realize that the universal regime sets in when the iner- 
tial driving of the collapse becomes dominant over the 
external driving force. This can be expressed by a local 
balance between inertia and the respective driving force. 
For system (i) the driving force is hydrostatic pressure 
and the relevant parameter thus the local Froude num- 
ber Fri oca i = r^/(gz c ) with gravity g and z c the depth 
below the surface where the cavity eventually closes. For 
system (ii) the local Weber number Wei oca i = pf^ro/a 
(with density p and surface tension a of water) gives the 
balance between inertia and surface tension as the rele- 
vant driving force. The duration of the universal regime 
can then be estimated as the time before pinch-off when 
these local quantities become of order unity. Figure [3J 
shows the local Froude and Weber numbers as a function 
of time to pinch-off r for a number of representative cases 
of system (i) and (ii), respectively. One can clearly ap- 
preciate that Weiocai for the needle system becomes unity 
later than Fri oca i for the impacting disc. This explains 
the large discrepancy in t u for the two systems. 

At the same time the distance between the two disc 
impacts with Fr=3.4 and Fr=4000 is much smaller than 
the distance between the two needle setups A and B. Ac- 
cordingly, the duration of the universal regime varies only 
between ~lms and ~6ms for the disc while its duration 
in the needle setup depends much stronger on initial con- 



ditions varying from micro- down to several nanoseconds 
as seen above. 

Our results allow us to explain the Froude-dcpendcncc 
of the experimentally and numerically observed exponent 
in [f| for the impacting disc. Based on Fig. [2] (c) these 
exponents can be viewed as a time-average of the local 
exponent. Due to the limited resolution and the onset of 
other effects (e. g. air flow) only the right part of these 
plots is accessible in experiments and the time-average 
will be heavily weighted towards the beginning of the 
universal regime. We can thus assume the experimen- 
tally observed effective exponent to be rougly equal to 
the maximum value of the local exponent. Using Eq. @ 
we can predict this value once the characteristic initial 
aspect ratio 7^ for each cavity is known as a function of 
the impact Froude number. 

This characteristic aspect ratio can be estimated to be 
proportional to the ratio between the (non-dimensional) 
maximum cavity expansion i? max and the closure depth 
z c . The first quantity scales with the Froude number ap- 
proximately as i? max ~ Fr 1 / 4 [l4| while the latter behaves 
as z c ~ Fr 1 / 2 [l4l. [24|. The characteristic aspect ratio of 
the cavity thus becomes 7; = R max /z c = CFr -1 / 4 with 
C a constant of order unity. Inserting ji into Eq. @ 
and solving for a (Ft), again neglecting da/dhir, we can 
thus predict the experimentally observable averaged ex- 
ponent which is found in excellent agreement with [f| as 
demonstrated by Fig. [3] (b). 

In conclusion, we have demonstrated that the universal 
theory of [ID, [23[ faithfully predicts the approach of the 
neck radius for invisicid, axisymmctric bubble pinch-off 
in four different systems widely studied in the literature 
over the past years. Remarkably, however, the duration 
of the final regime is shown to be strongly dependent on 
the type of system and the various control parameters 
employed. While it lies easily within experimentally ac- 
cessible time scales (~ms) for an impacting circular disc, 
it can be as low as a few nanoseconds for gas bubbles 
injected through a small orifice into a quiescent liquid 
pool. We were able to trace this difference back to the 
relative importance of the respective driving forces. Our 
findings reconcile the prediction of universality in bubble 
pinch-off [l2l . [l3| with an apparent dependence on ini- 
tial conditions [5||, an apparently constant scaling expo- 
nent @, 0, , and with the observation that non-inertial 
forces can be dominant in many experimental settings 
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